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Abstract. Here we treat the transmission of disease through a popu- 
lation as a standard Gahon- Watson branching process, modified to take 
the presence of vaccination into account. Vaccination reduces the num- 
^ , ber of secondary infections produced per infected individual. We show 

f^*) I that introducing vaccination in a population therefore reduces the ex- 

^P I pected time to extinction of the infection. We also prove results relating 

■^ , 

(^ the distribution of number of secondery infections with and without vac- 

o 

^ ■ Key words: Branching process, epidemic, vaccination, uniform conver- 

ince, gamma fi 
MSC: 60J85. 



X 



^ ' gence, gamma function. 

> 



1. Introduction AND BACKGROUND 

After the introduction of penicillin and mass vaccination campaigns of 
the 1950s, it was thought that infectious diseases would soon be a thing of 
the past. However, infectious diseases have turned out to be an intractable 
problem, as new infectious diseases continually emerge while previously 
existing diseases circumvent existing control measures through evolution. 
Accordingly, there is a continuing need to develop mathematically rigor- 
ously methods that can be applied to infectious disease epidemiology, in 
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2 ARNI S.R. SRINIVASA RAO* AND CHRIS T. BAUCH 

order to predict future spread of infections and assess the likely impacts of 
alternative control measures, including vaccination. 

One such area of promise emerging from the field of probability is branch- 
ing processes. Infectious disease transmission can be thought of as a branch- 
ing process in which each infected individual in a given generation of infec- 
tion produces some number of infected individuals in the next generation, 
as described by a probability density function. This picture is particularly 
apt in the early stages of an outbreak, when the proportion of infected indi- 
viduals is small and the process is highly stochastic. Epidemiologists need 
to know whether a disease will die out, and how many individuals are likely 
to be infected before the extinction, given certain fundamental probabil- 
ity distributions such as the distribution of number of secondary infections 
produced per infected individual. They also need to know how these param- 
eters will change once control measures such as vaccination are introduced. 

Some work has started in applying branching process-type arguments 
successfully to the mathematical modelling of disease transmission [1]. 
However much work remains to be done in developing a mathematically 
rigorous methodology to understand impact of vaccines. Branching pro- 
cesses have been applied to areas in biology such as the study of offspring 
distributions and the extinction of family surnames. However, there are 
issues that are particular to infectious disease epidemiology, such as the 
impact of vaccination on the branching processes. These pecularities ne- 
cessitate an extension of branching process methodologies to the case of 
infectious disease transmission and vaccination. Here, infectious disease 
transmission and vaccination are treated as a classic Galton-Watson branch- 
ing process. 

Let Zn be the number of infected individuals in the n''' generation of dis- 
ease transmission, where n > 0. Since we consider the early stages of dis- 
ease transmission and the proportion of susceptibles is almost 1, infected 
individuals do not co-infect the same susceptible individuals. An 'infected' 
(resp. 'susceptible') here means an individual who is infected (resp. sus- 
ceptible) in the present generation. Every infected individual has the same 
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probability distribution of the number of individuals they infect, and these 
are distributed independently of one other. Assume each infected in the n'^ 
generation generates A^ new individuals in the {n+ l)'^^ generation. Here 
A'^ G Z+is a measurable function (a r.v.) whose infected population distribu- 
tion is {Qi}T=o i.e. Pr [A^ = /] = 0/, / > 0. 

Let A^^ represents the number of infected by the /'^ individual in n''^ 
generation. These infected are account for the size of the population in the 
(n + l)"" generation. Then Z„+i = n["^^^ +A^f +^^ + ... + A^"^^^- Assume 



that each A^^ is iid r.v.'s i.e. Pr 



A^,(") 



= Pr[A^ = /] =0,-. The size of 
the {n + \y^^ generation Z„+i depends on the number of infected individu- 
als in the n''^ generation, therefore Z„ satisfies the properties of a Markov 
chain. We assume that Zq = 1, guaranteeing that there is one infected in- 
dividual from which the epidemic can start. After one disease generation 
(one discrete time interval), this individual has infected a further number 
of individuals as given by the probability distribution{0,}^Q. By this time, 
the original infected individual has also recovered (or died) from the infec- 
tion. Thus the first generation constitutes Zi infecteds, which is the sum of 
the Zq random variables each with a probability 0,. The second generation 
constitutes Z2 infected and the process continues. We call this process an 
epidemic branching process (EBP). 

Let (3„ be the proportion of the infected individuals who infect n individ- 
uals. Hence, |3o is the proportion of infected individuals who do not infect 
anyone before they recover or die. If Po = 0, then Z„ > Zq for all n and the 
process is not interesting, so we assume that (3o > 0. Likewise if |3o = 1 
then Kn, the probability of extinction of the virus after n generations, will 
be clearly equal to for all n. Hence, we assume < (3o < 1 throughout. 
Clearly, 7r„ < 7r„+i for all n. If the infection is extinct in generation n, it will 
remain extinct for n + 1 and all following generations. Then, {nn}n=i 2 
is a monotone increasing sequence, so we let n = lim„^oo7r„. Here n is 
the probability that the infection will eventually die out in the population. 
Suppose F{x) = po + |3i.x; + (32.x;2 + ... then F(7r„) = po + Pi7r„ + p27r,^ + ... 
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(b) 




Figure 1.1. Probabilistic realization of branching process 
{Z„}^^Q without vaccination (a) and with vaccination (b) 



= Y. ^rT^n- Fi'^n) is the suHi of the probabilities of the r + 1 distinct ways in 
which the infection can become extinct. 

In the next section we introduce probability generating function I{x) for 
{0,} and prove results similar to the one discussed above onF{x) in broader 
way. We study uniform convergence properties and related consequences of 
I{x) for X e (— M,]R) and x E ( — 1, 1). The results on uniform convergence 
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ofI{x) provided in this work are not new because I{x) is a power series (in- 
troduced in the next section). We study extinction probabilities with respect 
to the mean number of infectious individuals generated by starting Zq = 1 . 
This mean number is defined as /'(I). In the mathematical epidemiology 
literature, /'(I) would be referred to as the basic reproductive number. The 
process {Z„} is called supercritical if /'(I) > 1 and the process is called 
sub critical if /'(I) < 1. We have not seen rigorous mathematical treatment 
of the branching process for the spread of virus in the population (except 
elementary branching process applications for understanding measles out- 
breaks lUl and for virus HIV spread within the host cell dynamics [l4|]). Ex- 
tinction probability results in the pre-vaccination (i.e. without vaccination) 
scenario are derived from branching process studies available on offspring 
distribution, surnames and ot her gene ral branching processes analysis (see 
for example, flllflflliflHlillllJand fundamental results on uni- 
form convergence can be seen in many classical analysis books). Hence 
some results in section 2 do not bring new mathematics, but rather new ap- 
plications. In section 3, we bring novelty through theoretical arguments on 
EBP when vaccination is introduced in the population. By vaccinating the 
individuals, the resultant probability density functions of infection densities 
shifts to the left and extinction of virus is quicker. We also explore some 
consequences of vaccination. 



2. Elementary results on extinction probability 

We saw Pr [Nr = i] = Qi and A'^^ are independent discrete random vari- 
ables. Suppose I2 and Jr are probability generating functions of Z2 and 
A^^ for r = 1 , 2, 3 ■ ■ ■ . Then we know 7i = /i because Zi = A^i , 7i .J2 = h- 
Suppose distribution of A^^ are same (say some J), then I2 = [J] where 
h = hi^) = T4=o^^[^2 = i]^' is the probability generating function of Z2. 
In general let /„ is probability generating function of Z„. Since Z2= Ni + 
N2 + --- +Nzi , it follows that Fr[Z2 = k+\ /Zi =k\=Fr [Z-=i Ni = k+\]. 
We can prove h = h {h ) and /„ = /„_ 1 (/i ) L12J . Let I{x) = ZT=o Pr W = i] x' 
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(|;c| < 1) be the probability generating function for {9/}. We prove conver- 
gence properties of /(;c) and hence bring elementary results on extinction of 
probability of the virus from the population. 

Theorem 2.1. Suppose I{x) converges for \x\ < M, then for a 5 > 0, I{x) 
converges uniformly on [— M + 5, M — 5] andl{x) is continuous and infinitely 
differentiable on any open interval (— ]R,]R). 

Proof It is enough to show existence of/' i.e. differentiability of/, then it 
follows continuity of/. This is because, 



lim,^p[I{x)-I{p)] = lim,^/-^^^^—^[x-p]=l\x).0 = 

X — p 

Hence / is continuous. I\x) = Y^i^i^'^^ i\x\ < K)- Also, as ? ^ oo , it 



follows that limsupy |0,|= limsupy ? |9,|. This is due to fact that / ^ oo , 
it follows that ^/i = I. Since limsup^|9/| as n ^ oo, it follows that ^9( 
is convergent. Therefore, Y^^i^' ^^^ L'9/Jc' have interval of convergence in 
( — M, M) . Since I{x) is a power series, the required results is straight forward 
and can be seen in several classical analysis books. Suppose the series 
'Ln=o'^nX" converges for \x\ < M, and define g{x) = Yn^o'^nX" {\x\ < M). 
Then g converges uniformly on [— R + £, M — £] , for every £ > 0. Here g is 
continuous and differentiable in (— M, M) and g'{x) = Y.n=i nanx"^^^ . 

The above theorem established that the probability generating function of 
the infected distribution is differentiable over the real line if \x\ < M holds. 
However conventionally, we define I{x) iox\x\ <\. D 

Theorem 2.2. For a given £, 5 > 0, let «e(L9,) = {x G M : |;c - 1^9,| < £} 
and'B^{\) = {xgM: |;c-l| < 5} anJ9, < 9,+i/or {? = 0, 1,2-- ■}. Then for 
every Se(E9,), there exists a 'B^{\) with the property that for all x G ®s(l), 
x^\, it follows that I (x) G ®e(L9()- 

Proof. Let 7] = £9,. Consider the sequence {71},=o i 2 • Here Tq = Qq, Ti = 
9o + 9i, r2 = 9o + 9i + 92 and so on. We assume 9o > 0. Otherwise, if 
9o = the branching process will terminate. Since 9o > 0, and from the 
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hypothesis we have, 0i > Go > and 02 > 9i and so on. Also, we know 
0,- < 1 because {0,} forms a complete probability distribution. Therefore 
{7^} is convergent. We have the very important result by Abel, which states 
that when ^0„ is convergent and asx-^ 1, it follows that lim£0„.x:" = £0„ 
(—1 < X < 1). Hence from this result it follows that x E !S§(1)^ I{x) E 



Remark 2.3. The elementary theorem 1241 and the arguments in the proof 
given can also be found for studying extinction of surnames of families, 
extinction of particular offspring, particle extinction in a volume of gas and 
other biological applications. Otherwise, this theorem is very well known, 
but we state it here in the language of infectious disease extinction in the 
population. This will help readers to compare the above elementary result 
with the new results in the next section when a vaccine that prevents disease 
transmission is introduced into the population. 

Theorem 2.4. Given I{x) for x E [0, 1] and size of infected population in 
zeroth generation is one. If the process is sub critical or critical then 
limPr [Z„ =0] = 1 as n ^ o° and if the process is supercritical then there 
exists a 71(7^ 1) G [0, 1] such that I{tz) = K and as n ^ oa it follows that 
limPr[Z„=0] =71. 

Proof. Consider the closed interval [0, 1]. We saw in theorem ITT] that I{x) 
is continuous in [0, 1] and also I'{x) = X^?0,x' > 0. Therefore, I'{x) is strictly 
increasing in [0, 1]. If the system is sub critical then /'(I) < 1. This means 
I'{x) < 1 for X G [0, 1). This implies, J^ I\x)dx < f^dx^l- I{x) < 1 - 
X =^ I{x) > X. I{1) = Y,Bi = I. This means the curve of I{x) never touches 
I{x) = X in [0, 1) and at 1 it will be 1. Therefore, 7(1) = 1 has to be a unique 
fixed point for I in [0, 1] .Now it is left for us to show that limPr [Z„ = 0] = 1 . 
We can easily prove that {Pr [Z„ = 0] } is convergent. Let n be the limit of 
this sequence, so that as n ^ 0°, it follows that limPr[Z„ = 0] = n. Also, 
lim/(Pr[Z„_i =0]) = I{n) and I{k) = n. But we saw above that, when 
^'(1) < 1, then 7(1) = 1 is unique fixed point. Hence limPr [Z„ = 0] = 1. 
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Since I'{x) = ^/9,vv:'^^>0 for x E [0, 1], I'{x) is strictly increasing, then 
there exists aa E (0, 1) such that for Z? G (a, 1) we will have 

(2.1) 1 < I'{b)<l'{l) 
=^ dx < / l'{x)dx 

Jb Jb 

(2.2) ^\-b < \~l'{b)^b>I{b) foxbe{a,\) 

These kind of above arguments can be found in basic setting of branch- 
ing process (see M, (chapter 12)], M (chapter 4)], lill (chapter 1)], yj, 
(chapter 6)], llSl (chapter 0)]). Equation |2. II is a situation of supercritical 
process. When process obeys such property then our aim is to show that 
there exists two fixed points 71 and 1 in [0, 1] such that tt 7^ 1. Every value 
of b between a and 1 (exclusive), the value of line y = xh greater than the 
value of the curve y = I{x) for ;c G (a, 1) on the XY— plane. Therefore, this 
situation is not conducive for us to find a 71 such that /(7r) = n. If we bring 
a relation I{x) — x = for some x < I then prove that x is unique, then we 
are ready to compute extinction probabilities. Consider I{x) — x. We have 
7(0) — = 00 > and from equation lZ^ /(Z?) < Z? for Z? G {a, I). This means 
lix) —x<Qfoxx — b. Thus, 



I{x) —X 



> forx = 
< forjc = b. 



The intermediate value theorem says if any continuous function on a 
given closed interval with values ranging from negative to positive, then this 
function will be zero for some value in between the same closed interval in 
the domain. Therefore I{x) —x = for some x G {0,b). Let I (71) = n for 
n G (0, b) . We are sure that in (a, 1) there will be no fixed point. Suppose w 
is the another fixed point in (0, 1 ) then either w G (0, 7r) or w G (71, 1 ) . Since 
w is another fixed point, /(w) — w = 0, 71 is also a fixed point, so I{ti) — 71 = 
and 7(1) — 1=0. 1{x) — x is differentiable and continuous in (0, 1), so 
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(2.3) 
and 



/(wi)— wi = forO<wi<w 

I{W2)—W2 = forw <W2 <Tl 



(2.4) 



/(wi)— wi = forTT < wi < w 

I{W2)—W2 = forw < W2 < 1 

These equations 12 . 3 1 & lOl to aether contradict the fact that I{x) in strictly 
increasing in (0, 1). Therefore w = n. This means limPr [Z„ = 0]=K. Hence 
we have two fixed points n and 1 in x G [0, 1] . D 

Remark 2.5. If the mean number of infections exceed unity, then the virus 
from the population will go extinct in a finite number of generations with 
probability n. Suppose Zq = N then each of these individual's process can 
be treated independently. Thus limPr [Z„ = 0] = n^ and limPr [Z„ = oo] = 
1 — n^ . limPr [Z„ = oo] — > when A^ -^ ooeven if tt >> 0. This means virus 
will go extinct even if population at the zeroth stage has a large number of 
infected individuals. 

Proposition 2.6. Let {Zn^i^} is the size of the infected population during 
post vaccine scenario. If{Zn} is convergent then {Z„(|)} is also convergent, 
but converse is not true for (|) > 0. 

Proof. Suppose {Z„} is convergent. Then for 5 > 0, \Z,n — Zn\ < 5 whenever 
m, n > A^ G N. We have, 

< M+A^2 + ---+iVz„=Z„+i 

Therefore \Zm,^ — '^n,^ \ < 5 whenever m,n>N EN. Hence {Z„,,^} is con- 
vergent, n 
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3. Spreading is restricted by vaccination 

Suppose the process {Z„} is controlled by introducing vaccination into 
the population, which protects susceptible individuals from infection and 
therefore reduces the number of secondary infections produced per infected 
individual. We model the probability density function of the number of sec- 
ondary infections per infected person as a gamma function, so that it covers 
a wide range of epidemiologically -plausible scenarios. Vaccination will 
shift the peak of the probability density function to the left. Our aim here 
is to estimate the corresponding change in the time to extinction brought 
about by vaccination. 

Proposition 3.1. Suppose the time to extinction without vaccination is la 

and P[N = Wa] = Ta [wa,Ta,(Xa] and suppose the time to extinction with 
vaccination is x^ and P[N = wt] = F^iw^,, T^, at,]. If mean ofTy is less than 
the mean ofTa then x^ <Xa- 

Proof. Let w^, wl means of Ta and F^. We know that by vaccinating the 
population, the new number of infections generated by one infected is re- 
duced and hence the sizes of {Z,,} at each stage n is affected. We have 
Wh <Wb ^ Tt,CLb < TaOLa- Hence the time to extinct will be earlier with 
vaccine than that of without vaccination. D 

We discuss here some further consequences of vaccination. Suppose 
L",-o^«9w« = j'^ooWa d{Ta [waJa,<^a]) = TaUa- Here d{r[.]) is gamma 
distribution function. We have seen in the previous section that if £wa0v,;„ < 
1 then the chance of extinction is one and if Y.'^a^wa > 1 then the chance 
is n for 71 < 1. This also means limPr [Z„^ = 0] = 1 as n^ ^ oo. Suppose 
JZ,i,=o'^bQwb = J^^Wbd{rb[wb,Th,ab]) = TbUb- Please note Wa and Wb de- 
note infected individuals and G^a and 0w^ denote infected population dis- 
tributions for pre and post vaccination, {wj}, j = 1,2,3,... are individuals 
in post vaccination scenario. Here limPr[Z„^ = 0] = 1 as n^ ^ oo. Note 
that Z„^ < Z„^ and nb < na (by proposition 13.11) If we consider probabil- 
ity densities of these two scenarios on the XF— plane then the mean of the 
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infection distribution with vaccination shifts to left of the mean of the infec- 
tion distribution without vaccination i.e. 7i,a^ < Ta<Xa- Suppose {7}}, {cXj} 
for 7 = 1 , 2, . . . , n are set of parameters in the plausible range of post vaccine 
scenario. If T = Mm{ri,...,r„} anda = A/m{ai, ...,a„} then post vaccina- 
tion infection density will have one of these three parameters combinations: 
/. {T, a^} or II. {Tj, a} or ///. {T, a}. In the cases of / and II the infection 
density could be realistic, but III is not certainly an expected density. To 
avoid situation — ///, we choose a set Ej consisting of all possible combi- 
nations of parameters {Tj.aj} such that Y— density is never a decreasing 
function. 

The corresponding inequalities due to situations /, II and /// are Tttj < 
Tada, TjU < TaUa and Ta < TaOLa for some j such that post vaccination pa- 
rameter sets are in Ej. Since w^ < h^ it can be viewed that x^ < Xq. However 
the mean of Tj for any j is less than the mean of Fq. This implies, 

Pr[Z„^=0](asn„^oo) < Pr[Z„, = 0](as n^ ^ oo) 

Also, Pr[Z„^ = 0] < [Znj = 0] . Note that Wq, hj are not associated with the 
same time axis because pre and post vaccination situations cannot occur at 
the same time. In section 2, Z„,Z„+i,... are associated to the same time axis. 
Consider the sequence {Pr [Z„^ = 0] }^ ^q j j... • Note that the suffix a is used 
here to indicate that size of the population in pre vaccine scenario. We have 
7(0) = Pr [Zi = 0] = 00 > and also Go < 1. Then by the mathematical 
induction and monotonic property of / we can show that 

/(Pr[Z(„_i)^=0]) < /(Pr[Z„,=0])</(l) 
^Pr[Z„„=0] < Pr [Z(„+i)^ = 0] < 1 

Similarly by considering the sequence |Pr [Z„. = O] } _r. , 2 ^e can 
show that 
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7(Pr 



(n-\) 







< /(Pr[Z„^.=0])</(i; 



Pr[Z„^.=0] < Pr 



in+l) 



< 1 



Thus {Pr [Zfi^ = 0] } and {Pr [Z„ = O] } are convergent sequences. More- 
over {Pr [Z„ = O] } converges earlier due to vaccination. 

Theorem 3.2. Let Q„j be the probability that an individual infects Wj indi- 
viduals during post vaccination, then {Lj(Lw ^j9w;)}^ ^ {'LjTf){J^jCi?j)- 

Proof. Since {ri,...,r„} and {ai,...,a„} are non-zero real numbers, we 
have from the Cauchy-Schwartz inequality, 



;■ J j 

/oo 
w,-j(r,K-,r,-,a,])}2 

Thus this inequality realtes the mean number of infections and the pa- 
rameter values in post vaccination. Hence for modeling super critical and 
sub critical situations with respect to parameter values this inequality can 
be utilised. D 

Proposition 3.3. Suppose the mean distance between the peaks of the in- 
fection densities is d{Pj,Pa)- Here Pa{wa,Ta{wa)) are Pj{wj,Tj{wj)) are 

the peak points for the pre and post vaccination densities. Then T{W]) = 

1 

rM±[^{Ljd{Pj.Pa)y-{Wa-W]) 

Proof Let Pj{wj, Tj{wj)) be the point representing the peak of the density 
function after vaccination, in Ej. lfPa{wa, ^aiwa)) is the point representing 
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without vacc 
with vacc. 




2 3 4 5 6 7 

number of secondary infections 



Figure 3.1. Examples of infection peaks with (Pj) and 
without (Pa) vaccination. For the with-vaccination scenario, 
the average peak over many possible types is shown (see 
text). 



the peak of the density function before vaccination, then Euclidean dis- 
tance between these two points is d{Pj,Pa) (see Figure EU). The mean ( 
d{Pj,Pa)) of this function over all possible j is 



d{Pj.Pa) 



Ljd{Pj,Pa) 



For j possible types of post vaccine densities, let Pj{'Wj,r{'Wj)) be the 
point corresponding to the mean peak. O is a point which is at a distance of 
Wa — Wj from Pj to the right. Suppose d{Pa,Pi) is the distance from Pa to 
Pj. Then we know from Pythogoren right triangle principle that d{Pa,Pj) = 
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{OPaY + {OPjY- From general principles of means, we can easily verify 
the fact that d{Pa,P'j) = d{Pj,Pa). Therefore, 



^lY,d{Pj,Pa)\ = {Wa-Wj)^+{T{Wa)-T{W])y 



T{wj) = T{Wa) ± 



^\lld{Pj,Pa)\ -{Wa-Wj^^ 



This proposition relates the shape of the distribution of the number of 
secondary infections per infected individual before vaccination to the same 
after vaccination. This may be useful in situations where both distributions 
have small skew and Wa, Wj correspond to the expectation values of the 
distributions. D 



Corollary 3.4. Suppose Ia{^^) and Ij{o^) are variances of infections cor- 
responding to the distributions {0n^} and {e„.}. Then we have Ia{<3^) — 

l^{na-l',{l)fFr[N = na]-l^{nj-l'jil)fFr[N = nj] 

na nj 

= l^nlPr[N = na] + 2{{fj{l)f-{lUl)f}-l^njPr[N = nj] + 

"j 

m)f-m)y 



l,nlPr[N = n^]-l^njPr[N = nj]-2{{l'^{l)) + {l'j{l))}d^^^ + 

na Hj 

{(/^(l)) + (/;(l))}JI^I 

£n2pr[A^ = nJ-£n2pr[A^ = n,-]-{(/:(l)) + (/;(l))}j|^l 
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In the above /„(!) — I'Al) = J|a|. The above relation indicates that dif- 
ference between variability between infection process also depends upon 
the adjusted distance between the peaks described above. If the process is 
critical then the above difference is expressed as a straight combination of 
second central moments. 

Remark 3.5. 

£n«Pr[Z„^=0]<l ^ £n,Pr[Z„. = 0]<l 
=^ d{Pj,Pa) < v^ for some v < 1 

Again, 

^-1 |- 1 J E«;Pr[Z„^ = 0] > 1 (no impact of vaccination) 

1 E "a Pr[Z„^ = 0] < 1 (impact of vaccination) 

Since vaccine has a positive role in reducing the infection process such 
that as rij > N (for some large A^) it implies J^nyPr[Z„^. =0] < 1. When 
the process {Z„^} is supercritical then initially the process {Z„ .} could also 
be supercritical and eventually {Z„j} will attain sub criticality or criticality. 
Hence even if £ n^ Pr [Z„^ = 0] > 1 , eventually d> \. The larger the value of 
d the larger is the impact of vaccination and extinction occures much earlier 
with probability Tijiiij < Ti). 



4. Conclusions 

Branching process analyses have provided a rigorous perspective on nu- 
merous physical problems in the past. Here, we have introduced the foun- 
dations for application of branching processes to studying the early stages 
of epidemic spread in a population, both with and without vaccination. We 
modified the standard branching process to account for the effects of vac- 
cination, showing how vaccination will decrease the time to extinction. 
There remains much work to be done in epidemic branching processes. 
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particularly in the incorporation of other types of control measures (quar- 
antine) and in the incorporation of realistic population structures such as 
age-structure. 
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